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Abstract 

There has been found an exact solution of the mixed problem for 
Shrodinger's compact ?7(m)-vector nonlinear model with an arbitrary sign 
of the coupling constant. It is shown, that in case of m^3 there is a 
new class of solutions - mixed (/(cr+/i)-vector solitons with "inelastic" 
(changing the form without the energy loss) interaction at a>l and strict 
elastic - at (7=1. They correspond to the color complexes consisting of 
(T-bright and /i-dark solitons (cr+/i=m) and they can exist both in self- 
focusing and defocusing medias. The universal A'^-soliton formula for the 
attraction and repulsion cases has been obtained by the method of Hirota 
for the first time. 

PACS42.65.Tg, 42.65.Hw, 03.40.Kf 



0. The evolutionary system of coupled nonlinear Shrodinger equations 
(NLS-m) 

m 

iLjip* = ajk\ipk\'^i>j, j = T~m (1) 
fe=i 

ipj G C, Lj =8(^+1 Cjd^^; ajk,Cj G R 

arises within the weak coupling limits in various nonrelativistic models of 
the nonlinear field theory. The integrability conditions and the exact solu- 
tions of NLS are of a broad practical interest (nonlinear optics, plasma, 
ferromagnetism, hydrodynamics, Bose-Einstein atomic condensates etc.[l- 
6]) together with an academic one. A strict mathematical derivation of 
the two connected parabolic motion equations that are equivalent to the 
system (1) when m = 2 and ci — C2 is given in [7] where the self-infiuence 
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of different polarization waves in the nonlinear media of a tensor response 
has been studied for the first time. At (1) Cj parameters determine the 
dispersion, and matrix coefficients Ujk when j ^ k and j = k, determine 
nonhnear interaction and self-action of tpj fields correspondingly. Sub- 
ject to the implication of the variable and the sign of the parameters 
product sign{cjajk) = sc^ the system (1) at a classical level describes 
a spatial or time evolution of the m-component field in the nonlinear 
cube-medium[7,8]; at the quantum one - Bose-gas with m-color degrees 
of freedom [9,10] with attractive (aej, > 0) or repulsive {eek < 0) point 
interaction. 

The exact intograblc cases of the system (1), by Liouvill implication, 
are highly limited and require satisfaction of the rigid conditions in space 
of the controlled parameters Cj,ajk- On the bases of Zakharov's theo- 
rems about the additional motion invariants one may show (the proof 
will be given in a separate work), that when the conditions are met 
ttjk = ±akk,Cj = ±c the system (1) admits representation of a zero 
curvature (Lax-pair) and is embeddable into the scheme of the inverse 
scattering method (ISM) . The integrable reductions of NLS-m arising in 
this case form the family of solitons vector models with a unitary U{m) 
and pseudounitary U{m,n) group of symmetry. The known exact solu- 
tions of the given family of models, for example, those of the U (2)-vector 
model of Manakov [8] (Lo = = L2-, cr = 2, /U = 0) 

iLov; = + i^2i')v;, j = M (2) 

are "single-color" multisolitons: bright solitons {tpi,2 ~ secha) in a 
self-focusing medium (ae > 0) [8] and dark ones (V'1,2 ~ tanh/3) in a 
defocusing medium (ai < 0) [11,12]. Bright-vector [8] and dark-vector 
[11,12] soliton solutions that are constructed on the base of the tradi- 
tional boundary problems of the Z7(l)-scalar bright (■!/'(±oo) = 0) and 
dark (■!/'(±oo) = pexp(i0)) solitons [5,13] may be classed as trivial ones 
in the given sense. The exact solutions of the mixed vector solitons in the 
pseudo-Euclidean [7(1, l)-model [9] and in the Euclidean t/(l+l)-modcl 
(2) with a defocusing (ae < 0) nonlinearity have demonstrated [12] that 
they also interact (as ?7(l)-scalar solitons) in the elastic (trivial) fashion. 

For the first time it was shown in the given work that in the family of 
i7(m)-vector nonlinear models of Shrodinger (integrable reductions of the 
system (1)) 

iLoi)] = se^ iV'fcl^V'j*, j = (3) 
fc=i 

with the boundary conditions of the mixed density 



V'a(C,0 



-^P^e'^', (4) 



^0, Vm(C,0 

at the rn 3, there exist a class of exact solutions - "color" f/((T+/i)- 
vector solitons with a nontrivial interaction (intermodc interchange) . The 
conditions (4) imply that every degree of freedom ■i/)„(l < n ^ m) in the 
system (3) has its vacuum (condensate) of Z6rO Pn = or finite pn =^ 
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density with an asymptotic phase G,i. In this case m = cr + /i, and in 
the rest a and /i, take on the arbitrary values 1,2, ...,m. 

A''-soliton formula of the m-component system (3-4) depends definitely 
from the medium character {sign ae = ±1) and in this sense it is a universal 
one for both self-focusing (attractive, +1) and defocusing (repulsive, — 1) 
media of cube-nonlinearity. 

ISM in case of the system (3-4) faces the necessity of analysis of the 
(m + l)-sheet Riemann surfaces and thereupon the AT-soIiton solution is 
obtained by a more economical (in a mathematical sense) method of Hi- 
rota [14]. It is shown that an elastic (trivial) interaction of the mixed 
(bright and dark) (7(1 + l)-vector solitons [9,12] is a consequence of the 
A^-soliton solution factorability for the system (3-4) in a particular case 
of m = 2. In the general position situation the A'^-soliton solution is not 
factorizable and the interaction of the color multisolitons has the inelastic 
nature (changing the form at the energy conservation). There have been 
found special cases when the interchange between nonlinear modes does 
not arise and the A'^-soIiton scaterring turns out to be a factorizable one. 

1. Let's introduce the functions of Hirota 

Gj = HiPj, Gj €C,H €B., j = 

and make a transition {Lj Dj) from the linear operators -L to the 
bilinear D ones that are defined hereinafter as 

D {U- V)^{D U)V - U{D V). 

Hereinafter variables and ^ at (1), (3-4) will be attributed the sense 
of time t and coordinate x. 

Taking into account the scaling changes |ae|=2, Cj=c{>0), x—^x^/c, the 
system (1) (due to the concept of Hirota) forms a bilinear family of the 
compact ?7(m)-symmetry 

DiGj-H = 0, 3 = 1^ 

m 

D2H-H = 2SY.\ Gk\\ (5) 

k 

{Di = iDt + D2, 1)2 = bl- A). 

where A G R is an arbitrary parameter which will be defined lower; 

5 = signse. 

In the formalism of the bilinear operators for the functions Gj and H 
there is a representation in a series form due to the formal parameter e. 
Let's chose this representation in such a way that it was coordinated with 
the non-trivial boundary conditions (4), i.e.: 

oc 

G j = ^ £^''{g0tig2jjlJ. + £90ag2u+lSja) 
v=0 

oc 

H = ^ e^^te^; g^ = ho = gocr = 1, (6) 
Sag — Kronecker's symbol, j = 1, m. 
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It is obvious that functions Gj determine jn-component field (m=CT+^) 
in an arbitrary combination a and ^ (for example, a-bright solitons and 
/i-dark ones). We will get A''-soliton solutions of the bilinear system (5) 
following the standard scheme of Hirota {R ~ Gj, H; j = 1, m). 

s" s""^ 

R = Ro — y Ri — y R2 — ^ - - - — *■ Rn- (7) 
Let's lay the vacuum solution go/a = Pi^ exp(i0^), = ki^x — {k^+X)t 

m—a 

in a zero order in e. One cam determine A = —25 ^ from (5). In the 

physics of optical solitons the sign function 5 = signse defines a self- 
focusing {S = +1) or defocusing {S = —1) character of the nonlinear 
medium. We have the following in a zero order in e: 

N m — a 

= X] 7n ' exp(??n), Vn = (nX + i{Cn + 2<5 pI^^' 
n=l fx 

where 'jn'^ and Qn are arbitrary complex parameters. In a one-soliton 
(A'^ = 1) sector series (6) in the scheme (7) stop in the second order in 
£, in a two-soliton (A'' = 2) sector they stop in the fourth order in e etc. 
The solutions of the system (5) describing the propagation dynamics of 
A'^-solitons in two a and /u-sectors of U (m)-vector ■i/'-space (fixed on the 6 
order in e) is represented as: 



N N 
n=l ij 

N 

E anijlm'niVm+ ■ 
l,m 



ij Im 

N 

+e®E<im9r-^9^*+--- 



N N 

H=e°+ E ViVj [e^o<j + E 'ni'nm{s^aijim+ 

ij Im 

N 

-\-S ^ Clijlmqr'nQ'nr'^ • 
q,r 



(8) 



(9) 
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Here e = 1 (R.Hirota), fjn = exp(77n); 



CqiCg/CrjC-rm J q4-^i*-^l 

^ijlmqr — ~Z "I T T O'ijlm^qr \ \ 

SQjSsqm(sirSlr r*-*j*-^m 

fj, fj, fj, J. .7 /'in\ 

^...ij... ~ ^ij^---ij----> ^ij — ~~Z 7 Zjfj, = C,j — iKfj., [W) 

m — fi 

E 7f7j - . 

cr = l 

Om = + Cm, Cim = " Cm, O S O*, fc;^ £ R. 

It is quite evident from here that the two-soliton [N = 2) solution 
stops at the 4th order in e. At the same time one can assure that the 
formulas (8-9) correspond to the exact three-soliton {N = 3) solution of 
the system under discussion (5). The solutions of the higher order are not 
given here because of their awkwardness. 

2. One-soliton (A''=l) solution of Shrodinger's mixed J7(a+/^)-vector 
nonlinear model (3-4) from (8-10) takes the form: 



\i>M ) ^ 5<.M(l + <ie''^+''^) j 

where Hirota function is = 1 + ane^^"^^^, 

m—cr m—fj, 

ar/ = c?i(.^+ E pIM-')/ E h?\\ 

IX a 

o-ii = ^fiflii, -^11 = -cxp(2i0i^), 
01^ = arctan((Im(i — fcp)/Re(i). 



As one can see the i7(CT-|-/i)-vector soliton of the mixed color (11) is 
a dynamics-and-topological formation and in particular cases coincides 
with known earlier one-color bright-vector (-(/;{^}=0, 5=-\-l) [8] and dark- 
vector (i/){ct}=0, 5=—1) [11,12] solitons. However, there is a principally 
new thing, i.e. the fact what the exact solution (11), unlike the vector soli- 
tons [8,11,12] has its place in the system (3) both in the cases of attraction 
(self-focusing, 5=+l) and repulsion (defocusing, 5=—l). Besides that, one 
should mention that a universal color (7(o'+/i)-vcctor soliton (11) can be 
in several states determined by its dynamics-and-topological nature. It 
is convenient to interpret these states in the language of particles. Let's 
denote the admissible isotopic states of the color U (cr -|-/u)-vector soliton 
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by the symbol {a,jj,}, where a+n—m. Then by analogy with a quantum 
chromodynamics, the state with a mixed color {u^O, fi^O} may be sup- 
posed to be "aromatic" and the one with a mixing miss ({0, /i}, {a, 0}) - 
"no aromatic" (one-color). In the given analogy the J7(m)-vector soliton 
of a mixed color has the intrinsic structure and the existence of different 
states is natural for such a compound particle. For example, in case of 
Shrodinger's J7(5)-vector nonlinear model the solution (11) for one {3,2} 
of the 4 admissible ({1,4}, {2, 3}, {3, 2}, {4, 1}) aromatic states of the color 
U (3-|-2)-vector solitons has the form: 

V^i \ / Ai(tanhX-|-itani/>i)e^®i \ 
V'2 = yl2(tanhX + 'ttan02)e"''2 . (12) 

(V'3V'4V'5)' / V {BaB^BsYsechXe''^ J 

2 

= Pf COS (j)^,@i^ = ki^x - (kl -26 pl)t + 0^, 
4>fi = arctan((w — 2A;^)/u); 2X = u{x — vt — xo), 

s.=7^[(E4 + '5«V4)/E l7n']i 

H=l a=3 

2e = vx + (m^ - + 85 E Pl)t/2; 
M=l,2; a = 3,4,5; 

whore (/,=2Re("i and u=2Im(^i are soliton's reverse width and ve- 
locity. It is clear that a mixed L'"(5)-vector soliton (12) consists of two 
dark {ipi,'(p2) and three bright (V'3,'i/'4,'05) components (nonlinear modes) . 
It is also obvious that the number of all the [/(5)-vector soliton admis- 
sible states is equal to 6, however two of thorn are no aromatic (one- 
color) vector solitons: a bright-vector one {5, 0} in case of self-focusing 
medium (S = +1) and a dark-vector one {0,5} in case of defocusing 
medium (d = —1). The change of the medium refraction index An^ ~ 
+ |''/'2|'^ + • • • + I'/'sP induced by the tJji, . . . , ips components interac- 
tion may be calculated by the forward substitution of an explicit solution 
(12). However, there comes a universal formula from the bilinear system of 
equations (5) for the whole family of Shrodinger's U (m)-vector nonlinear 
models, 

An'=Y:pl + S-^lnH, (13) 

which permits to determine the value of An^ of the Hirota unified function . 
So, for example, in case of consideration of the above [7(5)-model fj, = 1,2, 
H = l+ oiie''i+*'i . There will be derived from the (13) 

An^ = Pi+ pl± ^-^^ sec/i^ 
for the self (de) focusing (-|-(— )) medium. 

One should mention that presence of vacuum-condensate in the system 

with a finite density in the defocusing (5=— 1) media imposes a natural 
limitation on the color vector soliton characteristics: "U^ ^ 4 ^ cos^ (j)fj,. 



u{x — vt — xa) 
2 
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Nevertheless as far as the number of the admissible aromatic states of the 
color [/(m)-vector soliton is equal to (m — 1), the detection of exactly 
the same states in the multimode optical systems may turn out an event, 
which is more probable than one-color states whose number is equal to 
two. 



3. Two-soliton {N = 2) solution and dynamics of the mixed color 
multisoliton interaction. 

Let's show that the interaction of the color multisolitons (8,9) in 
Shrodinger's mixed ?7(m)-vector nonlinear model (3,4) is nontrivial (chang- 
ing a form without any energy loss) at m ^ 3 and there has a place 
the intermode interchange (energy) that is proportional to the nonlinear 
modes intensity of solitons. Let's study for this goal (without community 
limitation) asymptotic {t — » ±oo) behavior of the color multisolitons (8,9) 
at Af = 2. 

Two-soliton {N = 2) solution from (8,9) takes the form (cr-|-//=m): 



(14) 



3 

^ 90^(1+ E aJ*^e''*+^J +a*'i22e''l+'5l+''2+*)2) 

where Hirota function 

2 

H = (l + J2 aije'''+^^ + aii22e''i+*^i+''^+*^^) 

i,j 

Let vi > W2(Im(^i>Im("2), where Vn is S'"-soliton velocity in j mode 
{j=l,2, . . . , n, . . . , n+a). Solution (14) at t — > ±oo on the paths 

= X — Vnt of separate solitons fall into the sum of the free one-soliton 
solution of the type (12): 



1p{j}ix,t) 



t — >zboo 

Here CJ^-amplitude, -^-envelope of j mode of n soliton n — 1, 2; j 



-- J2 C;*5{" } i^-vnt, xt)e'^"' (15) 



a, n: 

S"^y = sechY^, S"^y = tanhy„=^ + itan0„j^, 
Yn = tin (a; - v„t - Xo„)/2, 
</'nM=arctan((t;„ — 2fc,i)/w„), M„=2ReCn, t)n=2Im^„. 

Amplitudes C7"* of 5" solitons before (-) and after (-I-) the interaction 

are connected by the relations C"^ = S"C"^ , where S-special matrix, 
converting asymptotic at t —oo in the asymptotic at t — > -|-oo: 

Sa =Cl2(l-Sl72l)(l-SlS2)-i 

2Cl- = 7f (aii)-i5; = e'(2*=''--\ 
Clr =p^cos0i^e*('^i^'+"\ 

Si = C2l(l - S27r2)-'(1 - S1S2)5, (16) 

2Cr = af,,(aii22aii)-i,Sl = g-'f^'AiM-), 
Cf,- =p^cos02^e'(2'^i^+'^2^>; 
= f , SI = If, .2 = fe, IC12I = IC21I = 1. 
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Hence it is quite clear that solitons velocities 7j,i-motion invariants, and 
phases X^^ and amplitudes C"^ are not such. As 1 5^1 7^ 1 and ISJ^I = 1 in 
the situation of a general position, it is quite obvious that the interaction 
of the mixed color U{a+fj,)-YectoT solitons (14) has a nontrivial (inelastic) 
character. As a result of such an interaction between nonlinear modes 
there arises the intensity interchange ~ \S^\^. The intcrmode exchange 
initiates energy redistribution in the components (nonlinear modes) of 
the color vector solitons. However, the given interchange phenomena in 
the a- and fi- modes have their peculiarities: (j-modes interchange by the 
finite energy ~ \S^\'^ and maintain the sign; fi- modes maintain energy 
(IS^P = 1), but change their polarity and attain an additional phase jump 
~ (20nM + tt) as an interaction result. It is clear that jj,- modes interact 
only with different values of phases. 

Asymptotic analysis shows in general that the exchange between the 
components of a separate color soliton is not arbitrary (chaotic) but is cor- 
related with the commensurable changes in the components of all the other 
solitons. A nontrivial interaction between the color multisolitons (14) and 
the admissible scenarios of the intermode exchange in the system (3,4) are 

m m 

regulated by the laws of conservation: a) ^ \C"~f = ^ IC""*"!^ -total in- 

i i 

tensity of a separate soliton and 6) E(E|Cn^) = E(E IC]'"^!^) 

n j n j 

- full intensity of the all the solitons before (-) and after (+) interac- 
tion. Validity of the given laws may be easily seen from the asymp- 
totic formulas (16). Besides this, shifts of their inertia centers AX„ = 
X+„ - Xq- , AXn = {-ir+^2Q^ Inx, where 



(«+ E pJ/I^ImI )(''+ E pI/\^2.\^) E 7f7j7r7j 
g _ f jvr 

IS+^T," pJ/^iM«2M P "e " \t? P l7j P 

arising as a result of solitons interaction obey the condition of Sudzuki- 
Zakharov-Shabat (the law of solitons center of inertia conservation): 
CiiAXi + (^22^X2 = 0. The last one is a consequence of the value conser- 
vation Itot = / E \'>Pn\'^dx in time t. The laws of conservation and strict 

n 

formulas that are mentioned above permit to determine the interchange 
"kinetics" and admissible scenarios of the intermode switches in the sys- 
tem of the color multisolitons (14) by the quantitative fashion. However, 
let's pay attention to the factor e in (17). 

Due to the apparent many particle effects in e the shifts of the inertia 
centers of the color solitons AX„ don't accord with a conventional fear 
ture of factorization that is traditional for the ordinary solitons. Thus, 
e (in a concentrated form) points at a complex nature of interaction of 
the mixed U{a + /i)-vector solitons (14). In the situation of a general 
position (solitons parameters 7^,Cn, (%>), vacuum densities and com- 
ponent number a+fj,=m - arbitrary), the shifts AX„ of solitons because 
of £ can't be presented in a two-particle form, iV-soIiton scattering doesn't 
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come to the pair one and the interaction of mixed U(a+fi) -vector solitons 
at (T^2 is nontrivial (changing the form at the energy conservation). In a 
particular case of the hnear dependence of the parameters fi —^i jj =0 
their influence on e disappear, the vacuum contribution are balanced and 
the center shifts AJf„ admit two-particle representation. Consequently, 
the A'"-soliton solution is factorized and the soliton interaction becomes 
elastic (|5']'|=1, j=cr, /u). Besides this, from the (17) directly results that 
in a special case of Manakov's mixed U (2)-model [12] a=jj,=l the inter- 
action of the mixed ?7(l-|-l)-vector solitons is strictly elastic. In all the 
other cases at m>3 color L''(cr+/i)-vector solitons {m—a+jj.) interact by 
the nontrivial fashion and there is an energetic interchange between their 
nonlinear modes that is consistent with the above-mentioned laws of con- 
servation. 

In conclusion the authors express their gratitude to V.E.Zakharov for 
his attention to their work, S.V.Manakov and V. G.Marikhin for their 
useful remarks. The work was partially fulfilled at Landau ITP. 
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